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In an interesting recent combined experimental and 
theoretical paper by Wang et al. on strontium clusters, 
the experimental magic numbers at N—3A and 61 were 
assigned to the two structures shown in Fig.QJ 1 An ad- 
ditional strong magic number at N=82 remained unas- 
signed. As both the assigned Sr34 and Srgi structures 
are polytetrahedral (that is the whole space occupied by 
the cluster can naturally be divided into tetrahedra with 
atoms at the vertices) it seemed to us that a natural 
interpretation of the magic number at N=82 might be 
in terms of the third structure in Fig. ^ which was re- 
cently suggested as a particularly stable polytetrahedral 
cluster £ 

The assignment of the Sr34 and Sr@i structures was 
based on calculations using a model potential that Wang 
et al. parameterized.- However, the assigned structures 
were not lowest in energy, but calculations of the free 
energies of these structures (within the harmonic ap- 
proximation) indicated that they became most stable at 
T=129 and 220K, respectively. 

To test our hypothesis for Srs2 we decided to look at 
the structure and thermodynamics of this cluster using 
the same potential as in Ref. QJ. As for Sr34 and Sr@i the 
global minimum of this cluster is not polytetrahedral, but 
based on a Mackay icosahedron. However, using paral- 
lel tempering, a method which has proved particularly 
useful for probing low-temperature solid-solid transitions 
in clusters where achieving ergodicity can be particularly 
challenging, we found no evidence that the 82-atom poly- 
tetrahedral structure depicted in Fig. ^ ever became low- 
est in free energy. 

This failure led us to check the thermodynamics of Sr34 
and Sr6i presented in Ref. Q using parallel tempering. 
The resulting heat capacities (Fig. [2J are very different 
from Ref. QJ showing no evidence of large peaks in the 
heat capacity at the expected transition temperatures. 




2.6 
2.4 - 
2.2 
2.0 
^ 1.8 - 
^1.6 
1.4 
1.2 - 
1.0 



61 I 



34 



150 200 
temperature / K 



FIG. 1: Possible structures for the magic number clusters 
Sr34, Sr6i and Srs2 along with their respective point groups. 



FIG. 2: Classical heat capacities for Sr34 and Sr6i calculated 
from parallel tempering simulations. 



Furthermore, quenches from configurations generated by 
the simulations provided definitive evidence that there is 
no structural transition to the polytetrahedral structures 
of Fig. ^ ( or t° defective structures based upon them). 
For example, for Srgi at T=300K the Td minimum was 
found only once in 4000 quenches, compared to 59 times 
for the global minimum; even above the melting tem- 
perature of this cluster, there is no sign of the relative 
stabilities of the two structures reversing. We confirmed 
that these results do not reflect ergodicity problems by 
running simulations starting from the D§h and Td struc- 
tures: the clusters quickly escaped from these structures, 
seldom to return. 

To identify the cause of the differences between our re- 
sults and those of Ref. we repeated Wang et al.'s cal- 
culations of the partition functions of the global minima 
and the D^h and Td structures using the harmonic ap- 
proximation. However, we found the same scenario; i.e. 
structural transitions driven by the much greater vibra- 
tional entropy of the polytetrahedral structures. That 
the D§h and Td structures are observed much less in 
the simulations than one would expect from these cal- 
culations, therefore points to a breakdown of the har- 
monic approximation as the cause of the huge overes- 
timation of the partition functions of the D§h and Td 
structures. This unusual behaviour — the harmonic ap- 
proximation usually gives a correct qualitative picture of 
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the thermodynamics — is worthy of further investigation. 
It may be that the small mean vibrational frequencies of 
the D§h and structures are indicative of low barriers 
and basis of attractions that are small in extent. 

Our results do not necessarily imply that the structures 
in Fig. ^ incorrect assignments for the experimental 



magic numbers. Although our calculations with Wang et 
aZ.'s strontium potential are unable to provide any evi- 
dence in support of these assignments, this may instead 
reflect deficiencies in this potential and the difficulty of 
reliably modelling the structure of metal clusters. 
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